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ffi ABSTRACT 

The majority of the extragalactic sources yet detected at TeV photon energies 
belong to the class of "high frequency peaked BL Lacs" (HBLs) that exhibit a 
spectral energy distribution with a lower peak in the X-ray band. Such spectra are 
well described in terms of a log-parabolic shape with a considerable curvature, 
and widely interpreted as synchrotron emission from ultrarelativistic electrons 

■^j- outflowing in a relativistic jet; these are expected to radiate also in 7 rays by 

the inverse Compton process. Recently we have compared the X-ray spectral 
parameter distributions of TeV detected HBLs (TBLs) with those undetected 
(UBLs), and found that the distributions of the peak energies E p are similarly 
symmetric around a value of a few keVs for both subclasses, while the X-ray 
spectra are broader for TBLs than for UBLs. Here we propose an acceleration 
scenario to interpret both the E p and the spectral curvature distributions in 

/^ terms of a coherent and a stochastic acceleration mechanisms, respectively. We 

show how the curvature parameter 6~ 0.3 — 0.7 of the synchrotron X rays, that 
depends only on the latter acceleration component, can be related to the inverse 
Compton luminosity in 7 rays, so introducing a link between the X-ray and the 
TeV observations of HBLs. 

Subject headings: acceleration of particles - BL Lacertae objects: general - galax- 
ies: active - radiation mechanisms: non-thermal - X-rays: galaxies 



Introduction 



The BL Lac objects constitute a rare class of Active Galactic Nuclei (AGNs). Their 
observational features include: weak or absent emission lines, high radio and optical po- 
larization, superluminal motions, and a typical double-humped spectral energy distribution 
(SED, vF v ). Their continuum emission is dominated by non-thermal radiations from radio 
to 7-ray frequencies, that make them the most frequently detected class of extragalactic 
sources at TeV energies. The observed broadband emission is widely interpreted as arising 



in a jet of relativistic particles closely aligned to our line of sight (l.o.s.) (Blandford & Rees 



1978). 



In the widely accepted framework of leptonic Synchrotron Self-Compton (SSC) radi- 
ation, the low-energy bump is constituted by synchrotron emission from ultrarelativistic 
electrons accelerated in the jets; the high-energy component is due to inverse-Compton scat- 
tering of these synchrotron photons by the same electron population (e.g., Marscher & Gear 



1985 Inoue & Takahara 1996). 



The BL Lacs come in two subclasses: the "low-frequency peaked BL Lacs" (LBLs) in 
which the synchrotron peak falls in the IR-optical range, and the "high-frequency peaked 



BL Lacs" (HBLs) where it falls in the UV- X-ray bands (Padovani & Giommi 1995). To 
mark the HBLs detected at TeV energies from those undetected, we refer to the former as 



TBLs and to the latter as UBLs (Massaro et al. 2011a) 



A convenient description of the BL Lac SEDs has been suggested by Landau et al. 
(1986) in terms of a log-parabolic (LP) model, i.e., a curved, parabolic shape in a double-log 
plot. Subsequently, the LP model has been frequently adopted to describe the X-ray spectral 



continuum in HBLs (Tanihata et al. 2004 Massaro et al. 2004, 2008) as well as the TeV 



emission from the TBLs (Massaro et al. 2006 Aharonian et al. 2009); it has been also 



used Gonzalez-Nuevo et al. (2010); Giommi et al. (2011) to describe BL Lac SEDs in the 



sub-mm and infrared bands. 

Such curved spectra are known to arise both by synchrotron or inverse-Compton radi- 



ations from electron distributions featuring in turn a log-parabolic shape (Massaro et al. 



2004 Tramacere et al. 2007; Paggi et al. 2009). 



Recently, we carried out an extensive investigation of the X-ray synchrotron emission 
of both TBL and UBL subclasses, based on archival observations carried out by BeppoSAX, 



XMM-Newton and Swift between 1997 and 2010 (Massaro et al. 2011a). On adopting the 



LP model, the X-ray SED of HBLs is described in terms of 3 parameters: the SED peak 
position E p , its maximum flux S p evaluated at E p (or the corresponding peak luminosities 
L p ~ At[D 2 l S p in terms of the luminosity distance Dl), and the spectral curvature b around 



E p QTramacere et al. |[2007| |Massaro et al. ||2008[) 



the E p dis- 



Comparing the spectral properties of TBLs and UBLs, we found that: (i 
tributions of UBLs and TBLs are similarly symmetric around a value of few keVs for both 
subclasses; and (ii) the X-ray spectral curvature b is systematically lower in UBLs than in 



TBLs, so that the former feature narrower spectral shapes (Massaro et al. 2011a). 



In this Letter we compare the X-ray synchrotron luminosities L p of the TBLs and the 



UBLs, as derived from our previous analysis (Massaro et al. 2011a). Motivated by the 



observational results recalled above, we interpret both the TBLs and UBLs E p distributions 
in terms of a coherent electron acceleration scenario, and those of b as due to accompanying 
stochastic acceleration. Finally, we provide a relation between the X-ray spectral curvatures 
b and the IC luminosities. 

We use cgs units throughout this Letter and assume a flat cosmology with H Q = 72 km 



s" 1 Mpc -1 , Q 



M 



0.26 and ft, 



0.74 



,e.g. 



Dunkley et al. 2009). In the following, the 



parameters E p , S p and L p refer to the observer reference frame, while all unprimed quantities 
refer to the jet frame. 



2. Log-parabolic Synchrotron Spectra 

Log-parabolic electron energy distributions (PEDs; i.e., number of particles per unit vol- 
ume and Lorentz factor 7) are generally written in the form n (7) = Uq (7 / 7 P )~ 2_r ' log< - 7//7p \ 
where 7 is the electron Lorentz factor, n® the normalization, 7 P = (7 2 ) 1 / 2 the mean particle 
energy (i.e., the peak of 7 2 ^(7)) and r the electron curvature parameter. Such PEDs repre- 
sent the general solution of the energy and time dependent Fokker-Planck kinetic equation, 
that includes systematic (e.g., electrostatic or electrodynamic) and stochastic (e.g., turbu- 
lent) accelerations, together with radiative and adiabatic cooling as well as particle escape 



and injection terms (Kardashev 1962 Paggi et al. 2009 Tramacere et al. 2011). 



In general, both the peak energy E p and the peak luminosity L p of a synchrotron SED 
emitted by a curved electron distribution depend on 7 3p , the peak energy of the distribution 
7 3 n(j). But in the case of a log-parabolic spectrum 7 3p is proportional to 7 P itself to imply 



l3p ' Tp 



10 l/2r 



7 P 



1Q1/106 



(1) 



given that the relation between the curvature parameter r of the PED and the homologous 
b of the synchrotron SED reads simply b = r/5 (e.g., Massaro et al. 2006). Then for a 
typical TBL with b ~ 0.3 or UBL with b ~ 0.7 (Massaro et al. 2011a) the ratio of 7 3p to 7 P 
is always < 3. 



The synchrotron peak energy E p is generally proportional to 7JL B 5, where B is the 
mean magnetic field and S the beaming factor 5 = T~ l (1 — /3cos(?9)) _1 , with T the bulk 



Lorentz factor in the jet and i? its opening angle (e.g., Paggi et al. 2009 Tramacere et al. 



2009). In the case of a log-parabolic PED, this relation simplifies to the proportionality 



E p oc 7^ B 5 10 1 / r on account of Equation 1 (e.g. 



Massaro et al. 



the synchrotron peak luminosity L p scales proportionally to n{^p 



2010 



1%B 



2X4 



2011b). In addition, 



5 when considering 



variations of the PED, spectral shape (Tramacere et al. 2011). But, in the case of log- 



parabolic PEDs, for the product n(7 3p ) 7! = N e %f(r) obtains, where N e is the total 
number of emitting electrons, and f(r) is a spectral factor depending only on r and ranging 
from 2.5 to 1.3 for the typical values of spectral curvature of the TBLs and UBLs, respectively. 

We compare the L p distributions of the two HBL subclasses by performing a KS test, 
and find that they are indistinguishable at a confidence level of 99% (see Figure [l]). The 
L p distribution for the two HBL populations may still differ within one standard deviation. 
As complementary test, we simulate the two distributions of log(L p ) with the same number 
of events (i.e., 55 for TBLs and 76 for UBLs), adopting the statistical approach described 
in Massaro et al. (2011a). We considered two gaussian simulated distributions, having the 
same median (< log(L p ) >(UBLs) = 43.80 and < log(L p ) >(TBLs) 44.17), the same variance 
(cr l p (UBLs) = 0.06 and o-£ p (TBLs) = 0.11) of the observed distributions and spanning the 
same range of L p . Then, we measured the KS variable, -D K s,simui, between the two simulated 
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Fig. 1. — The L p distribution of UBLs (red) and TBLs (black). TBL data do not include 



Mrk 421 and PKS 2155-304 and the giant flares of Mrk 501 and 1H 1426+421 (Massaro et 



al. 2011a). The maximal separation D KS between the two cumulative distributions, (i.e, 



the KS test variable) is reported. 



distributions. We built the distribution of -DKS.simui repeating the simulations 8000 times and 
we found that the probability to obtain, randomly, the observed .Dks is 86%. Thus, UBLs 
and TBLs have the same E p and L p distributions, and differ, significantly, only as for the 
spectral curvature b. 

Accordingly, we can assume that j p , B, and 5 have close values in both TBLs and UBLs. 
Then, given the proportionality relations for L p and E p written above, these observational 
results allow to consider a similar number N e of emitting electrons for both HBL subclasses. 
Moreover, the impact of limited varia tions of L p is softened on 7 P and even more on B, since 



we expect B 



1/2 
7p 



see Section 



3.1 ); furthermore, the values of E p are closely the same for 

o in 

the two populations, and so L p turns out to be proportional to 7/ . 



3. Electron acceleration in BL Lac objects 

The characteristic energy 7 P of the PED, simply relates to the synchrotron SED peak 
energy by E p oc 7 p as said in Section pi is mainly set by the systematic acceleration com- 
ponent. On the other hand, the stochastic acceleration mechanism is responsible for the 
curvature r of the PED, related to the spectral width b. 

In the following, we assume the acceleration mechanisms effective in BL Lac jets to be 
a combination of systematic coherent acceleration, responsible for the energy peak position 
of the PED, and of stochastic acceleration, which accounts for the broadened PED around 
its peak related to the lower spectral curvature. 

We consider inverse Compton radiative losses to be subdominant compared to those by 
synchrotron emission, consistent with the observational lack of inverse Compton dominance 
in the HBL SEDs, (e.g., PKS 0548-322, [Aharonian et aF|2010| 1ES 0806+524, [Acciari et al. 



2008), for which any external component appears to be necessary to describe their spectral 



evolution. 



3.1. 



According to Cavaliere & D'Elia 



Systematic acceleration 

(2002), BL Lac jets are likely to be powered by 



the Blandford & Znajek (BZ, 1977) or the Blandford & Payne (1982) mechanism (see also 
Lovelace |1976 Ghosh fc Abramowicz |1997 Krolik |1999 Livio |1999 ), i.e., by the extraction 
of rotational power from a spinning supermassive black hole via the Poynting flux associated 
with the adjacent magnetosphere. 



As discussed by the above authors, the simple force-free condition E-B_ = governing 
these magnetospheres is likely to break down at the jet boundaries, due to considerable 
electric fields E < B parallel to magnetic field (e.g. Cavaliere & Morrison 1980); such fields 
are present in the BZ configuration, especially at the jet boundary. Alternatively, they may 



result from magnetic field reconnections in current layers at the jet boundary (e.g. Litvinenko 



1996, 1999); the related systematic acceleration mechanism is primarily electrostatic (see 



Massaro et al. 2011b for the case of gamma-ray bursts (GRBs). 



However, these electric fields will be electrodynamically screened out by the embedding 
plasma at distances that exceed the screening length d defined by 

2\ V 2 

<1 = — [ rp "": - 1 = 5.3 • 10 6 

J p 



% 



1/2 



cm 



(2) 



Air e 2 n J \ n 

( Cavaliere fc D , Eha~||2002 ). Here u p is the plasma frequency, j p is the characteristic electron 
Lorentz factor, m e is the electron mass, e its electric charge, c the speed of light, and n the 



total electron (e.g., Massaro et al. 2011b). 



Accordingly, the electron energy gained for each acceleration step writes 

7 a m e c 2 ~ e B d. (3) 

Assimilating 7 a = j p , from Equations [2] and [3] we obtain the expression for the typical Lorentz 
factor of the accelerated electrons 



B 1 



1p 



4 7r m e c 2 n 



9.8 • 10 4 



B 2 



?? 



(4) 



The peak energy E p of the synchrotron emission for an electron of Lorentz factor 7 p 
10 5 falls in the X-ray band, on adopting for the parameters standard values for HBLs 5 



25, n ~ 1 cm" 3 and B ~ 1 G (e.g., Massaro et al. 2006 Paggi et al. 2009 Celotti fc 



Ghisellini 2008; Acciari et al. 2008 Aharonian et al. 2009). 



The maximal energy radiated by an electron with Lorentz factor 10 5 is E ~ 0.05 TeV. 
Then, considering again 8 ~ 25 we expect the peak of the inverse Compton emission to lie at 
around E ~ 0.5 TeV. This is consistent with all TBL observations that have a 7-ray photon 
index typically 2 or higher in the TeV energy range, implying the energy peak of their inverse 
Compton component to lie below a few TeVs. 

To complete the above scenario, the limiting electron Lorentz factor 7^ attained by 
an electron corresponds to the condition where the acceleration compensates the radiative 
(mainly synchrotron) losses. This occurs for 



1M 



3e 

2n <r T B 



1/2 



1.9- 10 7 B- 1/2 , 



(5) 



at values considerably higher than given by Equation [4] for our standard values of B and 



n. 



Thus, the maximal energy available to the bulk of the electron population is of order 
7m Tn e c 2 ~ 1 TeV for B ~ 1G, in agreement with the observed TeV spectral "tail" of TBLs. 
Finally, we note that lower values of the magnetic field coupled with high Doppler factors 
have also been found to provide good fits of the HBL spectra during high luminous states 



e.g. Finkeet al. 2008). 



3.2. Geometry and timescales 

The simplest source condition obtains when the acceleration and the emitting region are 
cospatial. In particular, we suggest that the emission arises from thin sheaths of thickness 
A R that bound the jet with radius R 3> AR, as shown in Figure [2] There the particle 
density is low, and the screening length is sufficiently long as to allow the electrons to attain 
the required high energies to emit up to TeV energies as discussed above. 

In this scenario the relativistic aberration of light concentrates radiation isotropically 
emitted in the comoving frame into a cone with opening angle •& ~ T _1 <C 1, where V is the 
bulk Lorentz factor of the jet. On the other hand, only photons emitted within •§ around 
the l.o.s. will be detected by the observer. 

The typical delay time r^ between two photons emitted simultaneously in the comoving 
frame from different points on the jet surface is 

r (i = -[i-cos(r- 1 )]^2r 2 - (6) 

c c 

The variability timescale t v consistent with the delay time ta of the photons in the observer's 
frame implies an upper limit on the physical length of the emitting region of order 

I ~ 2Y 2 t v c . (7) 

This differs from the one usually adopted based on the motion of the jet bulk toward the 
observer and implying I ~ 8t v c. Such a "flashlight" effect is analogous to the one presented 
by Ryde & Petrosian (2002) in the case of GRBs, but adapted to a cylindrical geometry, 
and so more like a flash along the jet axis over I ~ 10 15 cm, produced by a relativistic shear 
instability. 

Observed variability timescales t' v of order 10 3 s imply I ~ R ~ 10 15 cm to hold for a bulk 
beaming factor 5 ~ 25. This is also consistent with the synchrotron loss length 

l s ~ 1.5 -10 19 j- 1 B~ 2 cm, (8) 



of order 10 14 cm for an electron with 7 p ~10 5 in a field B ~ 1 G. Considering these emitting 
regions to lie at the base of the jet close to the supermassive black hole central to the AGN, 
the jet radius R can match the height I (i.e., / < R, see Figure pi). Assuming for the thickness 
of the effective acceleration and emitting region AR ~10 -2 R, the total volume filled by the 
emitting electrons is V ~ 7r/2 ARIR ~10 44 cm 3 , close to the standard values considered for 



leptonic radiation processes of HBLs (e.g., Massaro et al. 2006 Aharonian et al. 2009). 
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Fig. 2. — A schematic view of the jet emitting region. 



3.3. Stochastic acceleration 



The observational evidence reported in Section [T] that UBLs feature systematically nar- 
rower spectra compared to TBLs may be interpreted in terms of a less efficient stochastic 
acceleration occurring in the former's jets. 

In fact, the curvature parameter of the PED is related to the stochastic acceleration 



term in a Fokker-Planck kinetic equation (e.g., Kardashev 1962 Tramacere et al. 2007 



2009 Paggi et al. 2009 Tramacere et al. 2011), and is inversely proportional to the 



stochastic acceleration rate; thus the synchrotron SEDs are relatively broader when the 
stochastic acceleration is more efficient. Specifically, the PED curvature r is directly related 
to the diffusion coefficient D in the Fokker-Planck kinetic equation by r oc D~ l . Higher 
values of D and faster diffusion also imply less time spent in the acceleration region (see 



Ginzburg & Syrovatskii 1964 Protheroe & Clay 2004). 



In a simple statistical picture, r is proportional to the energy gain e itself, while it is 
inversely proportional to the number n s of acceleration steps, and to the variance of of the 



energy gain; in sum, r oc e/(a^n s ) (Massaro et al. 2004; Tramacere et al. 2011). As TBLs 



and UBLs show similar E p distributions, we can assume that both subclasses have similar 
values of B and of the e/of ratio. Thus higher values of n s correspond to lower values of r; we 
suggest such high values of n s to be comparable with smaller acceleration regions, since each 
acceleration step is shorter. So, the observational evidence that b is systematically larger in 
UBLs than in TBLs is consistent with larger volumes V for the former than for the latter. 

Finally, we remark that while the average magnetic field is comparable in TBLs and 
UBLs, the difference is due to the small scale fluctuations in the power spectrum. On large 



scales B is related to the electric field (see Section 3.1 ) which is responsible for the systematic 



acceleration; on small scales, a turbulent component gives rise to stochastic diffusion (see 
e.g., Brunetti & Lazarian 2011, for a related approach concerning radio volumes). The latter 
component yields different numbers of acceleration steps in TBLs and UBLs, but averages 
out on large scales. 

On the other hand, considering that the similarities between the E p and L p distributions 
of TBLs and UBLs have been interpreted in terms of similar numbers of emitting electrons 
N e (see Sect. |2|, the curvature- volume relation described above suggests that the electron 
density n e = N e /V is larger in TBLs than in UBLs, making a brighter inverse Compton 
peak in a SSC scenario. We conclude that the TBLs, not only have wider X-ray spectra, 
but are also expected to be brighter in 7-rays than UBLs. 

In addition, the diffusion/acceleration timescale is inversely proportional to the diffusion 
coefficient; so low values of D will correspond to less variable sources. This feature is also 
consistent with the lower variability and the lack of giant X-ray flares found in the whole 



sample of UBLs in comparison with the TBLs (Massaro et al. 2011a). 



4. Conclusions and Discussion 



In Massaro et al. (2011a), we analyzed and compared the X-ray spectral properties of 



TBLs and UBLs, finding that they have similar E p distributions, both symmetric around a 
value of few keVs, while the X-ray spectral curvature b is systematically lower in the former. 

In this Letter, we have compared the L p distributions of the two HBL subclasses, finding 
them similar at high level of confidence level. Then UBLs and TBLs differ mainly as for the 
spectral curvature b; these observational results likely imply similar numbers N e of emitting 
electrons for both subclasses. 



We have proposed to interpret the E p and b distributions on assuming that the electron 



10 



energy gain is due to both coherent and stochastic particle accelerations. The scenario is 
based on re-acceleration rather than continuous injection of fresh highly relativistic electrons; 
re-acceleration occurs via both systematic and stochastic mechanisms, with equilibrium oc- 
curring between the overall acceleration rate and the radiative losses. 

Describing the coherent acceleration in terms of energy gain from an electric field, we 
have derived a relation for the expected particle Lorentz factors 7 P < 10 5 . Thus for a typical 
HBL with magnetic fields B ~1 G, plasma density n ~ 1 cm -3 and a beaming factor 5 ~ 
25, the expected synchrotron peak energy is at E p ~ 1 keV, as in fact observed in the X-ray 



SEDs of HBLs (Massaro et al. 2008, 2011a). 



On the other hand, the stochastic acceleration component is mainly responsible for 
spectral broadening around E p . In fact, the curvature b of the X-ray spectra is only depen- 
dent on the stochastic acceleration term in a Fokker-Planck equation, and thus is inversely 
proportional to the diffusion coefficient D and to the stochastic acceleration rate p acc , that 
is, b oc p~l c . Thus, we interpret the narrow X-ray SEDs of UBLs in terms of less efficient 
stochastic acceleration compared to TBLs. 

Finally, pursuing the stochastic acceleration scenario we have linked the curvature pa- 
rameter b to the volume of the emitting region, through its inverse proportionality to the 



number of acceleration steps n s (see Section 3.3). This curvature- volume relation, combined 



with the above consideration of similar values of N e , indicates the emitting electron density 
to be larger in TBLs than in UBLs, making the inverse Compton peak brighter in the former 
than in the latter. 

Thus electron energies sufficiently high to radiate in the TeV range are related to suffi- 
cient bright luminosities for effective detection. Conversely, narrower SEDs and lower fluxes 
make UBLs harder to detect in the TeV range than TBLs, in agreement with our previous 
results concerning their X-ray observations. 
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